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Abstract — A  mathematical  model  and  solution  procedure  for  predicting  the  thermal  performance  of  single¬ 
pass  solar  air  collectors  is  presented.  By  omitting  or  providing  a  top  glass  cover  over  a  plane  absorber 
plate,  or  by  providing  a  bottom  plate  under  the  absorber  plate  and  circulating  air  over  one  or  both  of 
the  air  channels  so  formed,  four  common  types  of  flat-plate  solar  air  collector  designs  were  considered. 
The  surface  temperatures  of  the  walls  surrounding  the  air  streams  were  assumed  uniform  whereas  the  air 
temperatures  were  assumed  to  vary  linearly  along  the  collector.  In  the  mathematical  model,  the  solar 
collector  was  assumed  sufficiently  short  for  which  the  assumptions  were  valid.  By  considering  a  steady 
state  heat  transfer  using  the  thermal  network  analysis  procedure,  a  set  of  simultaneous  equations  for  the 
mean  temperature  of  the  walls  and  the  air  streams  were  obtained.  Instead  of  solving  the  simultaneous 
equations  for  mean  temperatures  explicitly,  a  matrix  inversion  method  was  employed  using  a  standard 
sub-routine  programme.  Because  heat  transfer  coefficients  were  temperature  dependent,  a  set  of  mean 
temperatures  was  approximated  which  allowed  the  heat  transfer  coefficients  to  be  evaluated  as  a  first 
guess.  An  iterative  process  was  then  created  that  enabled  the  mean  temperatures  for  the  collector  to  be 
calculated.  The  newly-calculated  mean  temperatures  were  then  compared  with  the  initially-guessed  temper¬ 
atures.  The  iterative  procedure  was  repeated  until  consecutive  mean  temperature  values  differed  by  less 
than  0.01°C.  After  this,  another  section  of  collector  with  a  length  equal  to  the  previous  one  was  added  to 
the  end  of  the  first  collector.  The  temperature  conditions  at  the  inlet  of  the  second  section  were  assumed 
equal  to  the  outlet  temperature  conditions  of  the  previous  section.  The  iterative  procedure  to  determine 
the  mean  temperatures  was  repeated  for  the  next  section.  Additional  sections  were  added  until  the  required 
overall  length  of  collector  was  considered.  By  this  procedure,  predictions  of  mean  wall  and  air  streams 
temperatures  for  a  collector  of  any  length  could  be  obtained.  Although  only  four,  common  single-pass 
types  of  flat-plate  solar  collectors  are  considered  here,  the  solution  procedure  could  be  extended  to 
encompass  most  other  collector  designs. 


1.  INTRODUCTION 

A  solar  air  heater  system  consists  of  an 
arrray  of  interconnected  solar  heat  collectors. 
Most  systems  incorporate  a  once-through  or 
single-pass  type  of  forced  air  circulation  with 
cold  ambient  air  intake  and  solar-heated  hot 
air  discharge  to  the  working  space.  The  collec¬ 
tors  may  be  connected  in  series,  parallel,  or  a 
series-parallel  combination.  A  series-only  con¬ 
nection  would  result  in  a  long  but  higher 
temperature  system  whereas  a  parallel-only 
connection  would  result  in  a  short  but  lower 
temperature  system.  Apart  from  the  length 
effect,  the  performance  of  the  system  would 
also  depend  upon  the  air  circulation  rate,  physi¬ 
cal  design  of  the  collectors,  and  prevailing 
ambient  conditions.  Therefore,  a  knowledge  of 
the  operating  characteristics  of  the  collectors 
would  enable  the  system  performance  to  be 
optimised. 

There  are  numerous  works  on  both  experi¬ 
mental  and  theoretical  performance  of  solar  air 
heaters.  Recent  works  by  Parker  (1981), 
Vijeysundera  et  al.  (1982),  Than  and  Ong 
(1984),  Biondi  et  al.  (1988),  Duffie  and  Beckman 


(1991),  Verma  et  al.  (1992),  and  Parker  et  al. 
(1993)  considered  the  steady  state  heat  balance 
equations  and  involved  linking  plate  efficiency 
and  heat  removal  factors  with  air  flow  rates, 
surface  wind  heat  transfer  coefficients,  and  heat 
transfer  coefficients  between  the  moving  air 
streams  and  the  surfaces  forming  the  flow  chan¬ 
nels.  The  derivation  of  these  factors  are  quite 
complex  and  require  considerable  effort  especi¬ 
ally  when  dealing  with  the  numerous  design 
configurations  available  involving  different 
combinations  of  geometrical  details  and  air  flow 
patterns. 

2.  OBJECTIVE 

The  following  types  of  flat-plate  solar  air 
collectors  are  shown  in  Fig.  1:  Type  (I) — Single 
channel  design  with  single  air  flow  between  top 
glass  and  bottom  absorber  plate;  Type  (II) — 
Single  channel  design  with  single  air  flow 
between  absorber  and  bottom  plates  (no  glass); 
Type  (III) — Double  channel  design  with  single 
air  flow  between  absorber  and  bottom  plates; 
Type  (IV) — Double  channel  design  with  double 
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Fig.  1.  The  basic  solar  air  collector  types:  (a)  Type  I,  (b)  Type  II,  (c)  Type  III,  and  (d)  Type  IV. 


air  flows  between  top  glass  and  absorber  plate 
and  between  absorber  and  bottom  plates.  In  all 
the  designs,  the  absorber  plate  is  black-painted 
over  the  top  surface  and  the  bottom  of  the 
collector  could  be  insulated  or  left  bare. 

3.  ANALYSIS 
3.1.  Thermal  network 

The  thermal  networks  for  the  four  types  of 
solar  collectors  considered  are  shown  in  Fig.  2. 
For  Type  (IV),  the  following  heat  balance  equa¬ 
tions  are  obtained  from  the  thermal  network  at 
the  points: 

71:  Si  +  KM-TJ 

+  h%  (Tfi-T1)=[/t(T1-Ta)  (1) 
Tn:  h2(T2-Tn)  =  h1(Tn-T1)  +  Q1  (2) 
T2:  S2  =  h3(T2-Tt2)  +  h2(T2-Tn) 

+  ht23(T2-T3)  +  ht21{T2-T1)  (3) 
Tf2:  h3(T2-Tf2)  =  h4(Tf2-T3)  +  Q2  and 

(4) 


T3:  hA(Tt2  -  T3)  +  hr23(T2  -  T3)  =  Ub(T3  -  TJ 

(5) 

For  Types  (I)  and  (II),  eqns  (4)  and  (5)  are 
not  valid  and  the  heat  balance  eqn  (3)  is  replaced 
by  the  following  equation  at  the  point: 

T2:  S2  =  h2(T2  —  Tn )  +  ht21(T2  -  TJ 

+  Ub(T2-T&)  (6) 

For  Type  (II),  the  absorber  surface  is  non¬ 
transparent  and  hence  there  is  no  transfer  of 
solar  radiation  to  the  bottom  plate,  viz., 
S2  equals  zero.  In  this  case,  the  heat  balance 
eqn  (6)  is  replaced  by  the  following  equation  at 
the  point: 

T2:  0  =  h2(T2-Ttl)  +  hl21(T2-T1) 

+  Ub(T2  -  TJ  (7) 

For  Type  (III),  there  is  no  heat  extracted 
from  the  system  between  the  top  cover  and  the 
stagnant  air  cavity.  The  air  is  assumed  to  recir¬ 
culate  within  the  space  under  free  convection  in 
an  inclined  wall  cavity.  Hence  eqn  (2)  does  not 
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Fig.  2.  Thermal  network  for  solar  air  collectors:  (a)  Type  I,  (b)  Type  II,  (c)  Type  III,  and  (d)  Type  IV. 


exist  and  eqns  (1)  and  (3)  are  replaced  by  the 
following  at  the  points: 

71:  S.  +  KM-T^  +  h^-T,) 

=  [/t(71  -  TJ  and  (8) 
T2:  S2  =  h3(T2-Tt2)  +  hBC(T2-T1) 

+  hI23(T2-n)  +  hl21(T2-T1)  (9) 


3.2.  Temperature  distribution  along  flow 
direction 

Consider  the  heat  flows  into  and  out  of  an 
element  of  thickness  dy  in  a  single  air  stream  at 
a  distance  y  from  the  entrance  of  a  collector  as 
shown  in  Fig.  3.  The  temperatures  of  the  walls 
of  the  duct  forming  the  air  channel  are  assumed 
uniform  but  not  equal  to  each  other  for  each 
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Fig.  3.  Heat  balance  on  element  along  flow  direction. 


wall.  The  temperature  of  the  air  stream  at  a 
distance  y  from  the  collector  inlet  is  assumed  to 
be  at  T{.  The  temperature  is  assumed  uniform 
throughout  the  cross  section.  At  the  outlet  of 
the  section  (y  +  <5y),  the  mean  air  temperature 
is  equal  to  [T{-\-  (d7}/dy)  5y\.  The  heat  balance 
for  the  element  is  thus 

mCTt+QW5y  =  mC^Tt  +  ^SySj  (10) 

Simplifying  and  dividing  by  Sy  we  obtain 

d7}/dy  —  QW/(mC)  (11) 

By  assuming  that  the  useful  heat  transferred 
to  the  air  is  uniform  along  a  collector  of 
length  L,  the  mean  air  temperature  at  the  outlet 
can  be  found  by  integrating  the  above  equation 
to  give 

Tf0  —  Tfi  =  QWL/(mC)  (12) 

This  assumption  implies  that  the  air  temper¬ 
ature  varies  linearly  along  the  collector  as  shown 
in  Fig.  4.  From  observations,  this  assumption  is 
valid  for  “short”  collectors.  The  mean  air  tem¬ 
perature  is  then  equal  to  the  arithmetic  mean 

Tf  =  (Tti  +  Tfo)/2  (13) 


. 

U^J 

•  . 

Fig.  4.  Linear  fluid  temperature  variation  along  collector. 


The  useful  heat  transferred  to  the  moving  air 
stream  can  then  be  written  in  terms  of  the  mean 
fluid  and  inlet  temperatures  as 

Q  =  2mC(Tf -  TfJ/iWL)  (14) 

For  the  case  of  the  two-channel  flow  assumed, 
where  the  two  streams  are  at  the  mean  air 
temperatures  Tn  and  T{2 : 

Q\  =  2m1C(Tf)  —  Tn3)/(WL)  and  (15) 
Q2  =  2 m2C(Tl2  -  Tf2j)/(WL)  (16) 

By  introducing  Y1  =  2mlC/{WL),  and  T2  = 
2m2C/(WL)  into  the  above,  Ql  and  Q2  may  be 
expressed  as: 

Qi  =  TdTn-Ttu)  and  (17) 

e2  =  r2(Tf2-Tf2>i)  (18) 

3.3.  Mean  temperatures  matrices 
By  substituting  eqns  (17)  and  (18)  into 
eqns  (2)  and  (4)  and  rearranging  we  obtain  the 
following  expressions  for  the  Type  (IV)  collector 

(hi  +  hr2i  +  l7t)Ti  —  hi  Tfl  —  hl2i  T2  =  St  +  Ut  Ta 

(19) 

hi  Ti  -  (^  +  h2  +  Ti )Tn  +  h2T2=-Tx Til  l 

(20) 

—  ht2i  T\  —  h2Tfi  +  (h2  +  h3  +  ht2i  +  hl23)T2 

~h3Tt2-hl23T3  =  S2  (21) 
h3T2  —  (h3  +  /14  +  T2)Tf2  h^T3  =  —  T2 7}2ii  and 
(22) 

-  hl23  T2  -  h4  T(2  +  (hl23  +  Ub  +  h4)T3  =  Ub  Ta 
(23) 

For  Types  (I)  and  (II),  eqns  (22)  and  (23)  are 
not  valid  and  eqn  (21)  is  replaced  by 

-  K21  Ti  -  h2  Tn  +  (h2  +  hr2i  +  Ub)T2 

=  S2  +  C/bTa  (24) 

Further,  for  Type  (II),  the  term  S2  in  eqn  (24) 
equals  zero. 

For  Type  (III),  eqn  (20)  is  not  valid  and 
eqns  (19)  and  (21)  are  replaced  by 

(hnc  +  K21  +  Ut)Ti  -  (hnc  +  hin)T2  =  Si  +  Ut  Ta 
(25) 

-(hl2 1  +  hnc)Ti  +  (hnc  +  h3  +  ht21  +  K23)T2 

-h3T[2-ht23T3  =  S2  (26) 
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In  general,  the  above  eqns  (19)-(26)  may  be  displayed  in  a  5  x  5  matrix  form: 


m(hl  +  hr2i  +  Ut) 

-hi 

_  hr21 

0 

0 

~Ti~ 

hi 

-(ht  +  ha+rj 

h2 

0 

0 

Tn 

—Kn 

—h2 

(h2+h3+ht21+hr23) 

— h3 

—  hr23 

t2 

0 

0 

h3 

-(h3+h4+r2) 

h4 

Tn 

0 

0 

— hr23 

-K 

(hr23  +  Ub  +  h4)_ 

.t3_ 

~UtTa+S~ 

-r  tTfU 

=  s2 

-r2rf2ji 

.  UhTm  . 

For  Type  (III)  the  equivalent  matrix  is  4  x  4: 


■(h„c  +  hr21  +  C/t) 

(h„c  +  hr2i) 

0 

0 

'Ti' 

■c/tT.+sr 

—  (h„c  +  hr2i) 

(hnc  +  h3  +  hr21  +  hT  23) 

—  h3 

—  hr23 

t2 

s2 

0 

h3 

-(h3  +  h4  +  r2) 

h4 

Tf2 

-r2Tf2>i 

0 

—  hr23 

— h4 

(hr23  +  ^b  +  h4). 

,t3_ 

L/bTa 

For  Types  (I)  and  (II)  the  equivalent  3x3  matrices  are: 


"(hi  +  hr21  +  I/t)  fcr21 

Tj" 

'ux+s  1" 

hi  -(hi  +  h2  +  ri)  h2 

Tn 

= 

-r,rfu 

—  K21  —  h2  (h2  +  hr2i  +  I/b) 

_t2_ 

s2+c/bra 

Further,  for  Type  (II),  S2  equals  zero.  In 
general,  the  above  matrices  may  be  displayed  as 

[A][T]  =  [B]  (27) 

The  mean  temperature  vector  may  be  deter¬ 
mined  by  matrix  inversion  using  standard  com¬ 
puter  packages 

[T]  =  [A]-1[B]  (28) 

3.4.  Radiation  heat  transfer  coefficients  from  top 
and  bottom  surfaces 

The  radiation  heat  transfer  coefficient  from 
the  top  surface  to  the  sky  referred  to  the  ambient 
temperature  can  be  written  as 
hrs  =  oefTi  +  Ts)(Tf  +  TS2)(T  -  T8)/(T,  -  Ta) 
(29) 

The  sky  temperature  obtained  from  Swinbank 
(1963)  is  given  by 

Ts  =  0.0552Tt5  (30) 


3.5.  Convection  heat  transfer  from  cover  due  to 
wind 

The  convection  heat  transfer  coefficient  due 
to  wind  from  McAdams  (1954): 

fcw  =  5.7  +  3.8  V  (31) 

Watmuff  et  at.  (1977)  give  this  equation  as: 

hw  =  2.8  + 3.3  V  (32) 

The  wind  coefficient  could  also  be  determined 
from  the  correlation  developed  by  Sparrow 
etal.  (1979): 

Nuw  =  0.0158  Re0,8  (33) 

where  the  Reynolds  number  is  calculated  based 
on  the  equivalent  diameter  Dh  equal  to  4  x  plate 
area/plate  perimeter. 

Table  1  shows  values  of  wind  coefficients  cal¬ 
culated  from  the  correlations  presented  by 
McAdams  (1954),  Watmuff  et  al.  (1977),  and 
Sparrow  et  al.  (1979).  Generally,  the  values 
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Table  1.  Wind  coefficients,  6W  (Wm  2K  ') 
Wind  velocity  V  (m  s-1) 

0.0  0.5  1.0  2.0  3.0 

McAdams  (1954)  5.7  7.6  9.5  13.3  17.1 

Watmuffef  al.  (1977)  2.8  4.3  5.8  8.8  11.8 

Sparrow  et  al.  (1979)  0.0  3.0  4.2  6.0  7.3 


obtained  using  McAdams’  correlation  was 
roughly  twice  that  given  by  the  other  two  corre¬ 
lations.  Unless  otherwise  stated,  McAdam’s  cor¬ 
relation  will  be  employed  in  this  article. 

3.6.  Radiation  heat  transfer  between  parallel 
plates 

The  radiation  heat  transfer  coefficients 
between  two  infinite  parallel  plates  set  1-2  and 
2-3  are  given  as 

K21  =  °{T\  +  7l)(7i  +  T2)/(l/£l  +  l/e2  -  1) 
(34) 


Nu,*,  =0.018  Re0  8  Pr0  4  (38) 

for  9500  <  Re  <  22,000. 

For  situations  involving  a  large  property  vari¬ 
ation,  the  Sieder  and  Tate  (1936)  equation  may 
be  employed 

Nu  =  0.027  Re0,8  Pr1/3(^//rw)014  (39) 

for  flow  in  smooth  pipes  where  Re  >  10,000  and 
L/D  >  60. 

A  more  accurate  correlation  which  is  also 
applicable  for  rough  ducts  was  developed  by 
Petukhov  (1970) 

Nu  =  (Re  Pr/X)(f  /8)(p/py,)n  (40) 

where,  for  104  <  Re  <  5  x  106,  X  =  1.07  +  12.7 
(Pr2/3  — 1)(//8)1/2;  n=0.11  (heating  with  uni¬ 
form  7^),  n  =  0.25  (cooling  with  uniform  Tw), 
and  n  =  0  (uniform  wall  heat  flux  or  gases). 

The  friction  factor  for  smooth  tubes  for  the 
above  equation  may  be  obtained  from 

/  =  [1.821og(Re)  —  1.64]  ~2  (41) 


and 

Kn  =  <y(T\  +  ri)(r3  +  r2)/(i/e3  +  i/e4  - 1) 
(35) 

3.7.  Forced  convection  heat  transfer  coefficient 
between  parallel  plates 

3.7.1.  Laminar  flow  region  (Re  <  2300). 
Heaton  et  al.  (1964)  proposed  the  following 
empirical  correlation  for  the  local  Nusselt 
number  for  laminar  flow  between  two  parallel 
flat  plates  with  one  side  insulated  and  the  other 
subjected  to  a  constant  heat  flux: 


Nu  —  NUoo  + 


a  [Re  Pr(Ph/L)]m 
1  +  6  [Re  Pr(Dh/L)]" 


(36) 


where  the  constants  are  a  =  0.00190,  b  = 
0.00563,  m  —  1.71,  n  =  1.17,  and  Nu*  =  5.4  for 
Pr  =  0.7. 


3.7.2.  Transition  flow  region  (2300  <  Re  < 
6000).  Hausen  (1943)  presented  the  following 
empirical  correlation  for  the  average  Nusselt 
number  between  the  beginning  of  the  heated 
section  and  the  position  L  for  flow  in  a  tube: 

Nu  =  0.116(Re2/3  -  125)  Pr1/3 

x[l  +(Dh/L)2/3](/r//iw)014  (37) 

3.7.3.  Turbulent  flow  region  (Re  >  6000).  Tan 
and  Charters  (1970)  recommended  that  for  flat- 
plate  solar  air  heaters,  the  heat  transfer  in  the 
developed  flow  region  should  be  obtained  from 
the  relationship: 


3.7.4.  Entrance  region  (L/D  <  60).  In  the 
entrance  region,  greater  heat  transfer  coefficients 
may  be  expected.  In  their  studies  on  the  effect 
of  the  thermal  entrance  region  on  turbulent, 
forced  convective  heat  transfer  for  an  asymmet¬ 
rically  heated,  rectangular  duct  with  uniform 
heating,  Tan  and  Charters  (1969)  determined 
an  expression  for  the  mean  Nusselt  Number 

Nu  =  Nu00(l+SDh/L)  (42) 

where  the  factor  S  was  found  to  be  adequately 
represented  by 

3.57  <  L/Dh  <60:  S  =  14.3  log(L/Dh)  -  7.9 
(43) 

L/Dh>  60:  S  =  17.53  (44) 

It  should  be  noted  here  that  S  is  0  when  L/Dh 
equals  3.57.  The  range  of  Reynolds  Numbers 
covered  in  their  studies  was  from  9500  to  22,000. 
Nusselt  (1931)  recommended  the  use  of 

Nu  =  0.036  Re0  8  Pr1/3(D/L)0  055  (45) 

in  the  entrance  region  for  10  <  L/D  <  400. 

For  noncircular  ducts,  the  tube  diameter  is 
replaced  by  the  hydraulic  or  equivalent  diameter 
defined  as 

Dh  =  4AJP  (46) 

For  flows  in  between  parallel  flat  plates,  the 
equivalent  diameter  Dt  is  twice  the  separation 
distance  between  the  plates. 

Figure  5  shows  the  Nusselt  Numbers 
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Length/Equivalent  diameter  (L/D) 


Tan  &  Chart  (Re  =  10000) 

Nusselt  (Re  =  10000) 

^Tan  &  Chart  (Re=22000) 

Nusselt  (Re =40000) 

-fpethukov  (Re=10000) 

+  Pethukov  (Re =40000) 

A  Hausen  (Re=2100) 

A  Hausen  (Re=6000) 

X  Colburn  (Re= 10000) 

-X- Colburn  (Re=40000) 

▼  Dittus  &  B  (Re= 10000) 

▼  Dittus  &  B  (Re =40000) 

A  Sieder  &  T  (Re=  10000) 

A  Sieder  &  T(Re=40000) 

+  Nusselt  (Re=22000) 

Fig.  5.  Nusselt  Number  versus  L/D. 


obtained  from  the  above  empirical  correlations 
for  comparison.  A  discontinuity  in  the  curve 
can  be  observed  in  the  Tan  and  Charters  ( 1969) 
correlations  at  L/D  >  60  because  of  the  maxi¬ 
mum  value  of  S  imposed.  As  a  result,  this 
correlation  is  not  convenient  to  use  in  the 
present  work.  It  is  better  to  use  Nusselt’s  (1931) 
correlation  in  the  entry  region  because  the  equa¬ 
tion  takes  the  entry  length  into  account  without 
any  discontinuity  in  evaluating  the  Nusselt 
Number.  The  values  obtained  using  Nusselt’s 
correlation  are  higher  than  those  given  by  the 
other  correlations.  By  comparing  the  Nusselt 
numbers  in  the  fully  developed  region 
(L/D  >60)  for  turbulent  flow,  the  ratios  of 
Nusselt  Numbers  from  Tan  and  Charters/ 
Nusselt  and  from  Petukhov/Nusselt  correlations 
are  obtained  and  tabulated  in  Table  2. 


Differences  can  be  expected  in  the  Nusselt 
Number  depending  upon  which  correlation  is 
employed.  The  variation  may  range  from 
±50%.  In  the  present  work,  the  heat  transfer 
coefficient  correlations  as  listed  in  Table  3  will 
be  assumed.  These  correlations  will  be  employed 
pending  better  ones  in  order  to  demonstrate  the 
qualitative  accuracy  of  the  theoretical  model 
employed  to  simulate  the  thermal  performances 
of  the  different  types  of  solar  air  heaters. 


3.8.  Natural  convection  between  parallel  plates 
Hollands  et  al.  (1976)  gives  the  following 
relationship  between  the  Nusselt  Number  and 
Rayleigh  Number  for  free  convection  between 
inclined  planes  lying  in  the  range  0  <  Ra  <  105 


Table  2.  Ratio  of  Nusselt  Numbers 


Ratio  of  Nusselt  Number 


Re  =  10,000  Re  =  22,000  Re  =  40,000 


Tan  and  Charters  (1969):  Nusselt  (1931)  0.79  0.79 

Petukhov  (1970):  Nusselt  (1931) _ 075 _ 068 
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Table  3.  List  of  correlations  employed  for  forced  convection  HTC 


Flow  Regime 

Range  of  Re 

Equation 

Laminar  flow 

Re  <2,100 

Heaton  et  al.  (1964),  eqn  (36) 

Transition  flow 

2,100  <  Re  <10,000 

Hausen  (1943),  eqn  (37) 

Turbulent  flow 

10,000  <  Re 

Nusselt  (1931),  eqn  (45) 

and  0°  <  <j>  <  60°: 


is  defined  as 


Nunc=  1  +  1.44  1- 


[‘ 


1708  (sin  1.8  <t>)u 


Ra  cos  (f> 


"  1708  i+  r  / 

Ra  cos  $ J  |_\ 


Ra  cos  </> 
5830 


(47) 

The  equivalent  diameter  employed  in  the  above 
definition  for  the  Nusselt  Number  is  the  separa¬ 
tion  distance  between  the  plates.  Also,  the  nota¬ 
tion  [  ]  +  is  used  to  denote  that  if  the  quantity 
in  the  bracket  is  negative,  it  should  be  set  equal 
to  zero.  The  natural  convection  coefficient  is 
specified  for  Type  (III)  for  the  air  cavity  between 
the  glass  and  absorber  surfaces. 

3.9.  Overall  top  heat  transfer  coefficient 

The  overall  top  heat  loss  coefficient  may  be 
obtained  from 

Ut  =  hw  +  hrs  (48) 

3.10.  Conduction  heat  transfer  from  bottom 
surface 

The  bottom  heat  loss  coefficient  is  given  by 
Ub  = - - - - -  (49) 

<  =  I(xbi//Cbi)+  VK 

In  the  absence  of  bottom  insulation,  the  insula¬ 
tion  thicknesses  are  set  equal  to  zero. 


n,  =  miC(Tfho  -  TJ/(HWL)  (52) 
and  for  the  second  air  stream,  if  it  exists,  is 

r>2  —  m2C(Tt2o  —  Ta)/(HWL)  (53) 

In  the  above  equations,  the  inlet  air  temper¬ 
atures  are  assumed  to  be  at  ambient.  For  the 
Type  (IV)  collector,  the  combined  instantaneous 
heat  collection  efficiency  of  both  the  air  streams 
is  obtained  by  adding  the  individual  air  stream 
efficiencies 


n0  =  n1  +  n2  (54) 

The  theoretical  model  assumes  that  air  tem¬ 
peratures  vary  linearly  over  short  collectors. 
Instantaneous  efficiencies  are  expected  to 
decrease  as  the  length  of  the  collector  increases. 

3.13.  Physical  properties 
The  physical  properties  of  air  are  assumed  to 
vary  linearly  with  temperature  TC  because  of 
the  low  temperature  range  to  be  encountered. 
The  following  are  assumed  for: 

viscosity 

p  =  [1.983  +  0.00184(7- 27)]  1(T5  (55) 
density 

p  =  1.1774-0.00359(7-27)  (56) 

thermal  conductivity 

k  =  0.02624  +  0.0000758(7- 27)  (57) 


3.11.  Solar  radiation  specific  heat 

The  solar  radiation  heat  flux  absorbed  by  the  C  =  1.0057  +  0.000066(T  —  27)  (58) 

first  surface  is 


51  =  a1H  (50) 

and  the  solar  radiation  heat  flux  absorbed  by 
the  second  surface  is 

52  =  m2H  (51) 

For  Type  (II),  S2  equals  zero  because  t  for  a 
solid  absorber  surface  is  equal  to  zero. 

3.12.  Instantaneous  efficiency 

The  instantaneous  heat  collection  efficiency 
of  the  first  air  stream  for  a  collector  of  length  L 


4.  THEORETICAL  SOLUTION  PROCEDURE 

The  theoretical  model  assumes  that  for  a 
short  collector,  the  temperatures  of  the  walls 
surrounding  the  air  streams  are  uniform  and 
the  temperatures  of  the  air  streams  vary  linearly 
along  the  collector.  A  long  collector  can  be 
assumed  to  be  divided  equally  into  a  finite 
number  of  short  collectors,  or  sections.  The  wall 
and  mean  air  temperatures  of  the  first  section 
are  initially  guessed  and  specified.  In  all  cases, 
they  are  set  equal  to  the  ambient  temperature. 
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Heat  transfer  coefficients  are  evaluated  accord¬ 
ing  to  the  initially-guessed  temperature  values. 
An  iterative  process  is  then  created  and  the 
mean  temperatures  for  the  section  calculated 
using  the  equations  derived  by  employing  a 
standard  package  matrix-inversion  subroutine. 
The  newly-calculated  temperature  values  are 
then  compared  with  the  previously  assumed 
ones.  The  iterative  process  is  repeated  until  all 
consecutive  mean  temperatures  differ  by  less 
than  0.01  °C. 

Another  section  of  collector,  with  length  equal 
to  the  previous  section,  was  added  to  the  end 
of  the  first  section.  The  mean  wall  and  air 
temperatures  of  the  second  section  of  collector 
are  then  set  equal  to  the  mean  wall  and  air 
temperatures  of  the  section  before  it.  The  inlet 
air  temperature  of  the  second  section  is  set  equal 
to  the  outlet  air  temperature  of  the  first  section. 
The  iterative  procedure  is  repeated  until  all  the 
sections  of  the  given  collector  are  considered. 
By  this  procedure,  wall  and  mean  air  temper¬ 
atures  can  be  predicted  for  the  complete  length 
of  collector. 

The  computer  programme,  based  on 


FORTRAN,  is  outlined  in  the  Flowchart  shown 
in  Fig.  6.  The  start  of  the  programme  considers 
the  first  section.  An  initial  guess  of  the  mean 
wall  temperatures  and  mean  air  streams  temper¬ 
atures  are  made.  For  the  first  section,  these 
temperatures  are  specified  equal  to  ambient 
temperature.  (For  subsequent  sections,  the  ini¬ 
tial  temperatures  are  specified  equal  to  the 
respective  temperature  values  of  the  respective 
previous  sections.)  An  iterative  process  is  then 
initiated.  The  programme  calculates  all  the 
required  heat  transfer  and  heat  loss  coefficients 
based  on  the  initially-guessed  temperatures.  The 
matrices  [A],  [T],  and  [B]  are  then  set  up.  A 
standard  matrix  inversion  subroutine  is  then 
called  to  invert  [A]-1  to  determine  a  new  set 
of  temperature  matrix  [T'].  Each  new  temper¬ 
ature  value  in  the  matrix  [T']  is  then  compared 
with  the  corresponding  initially-guessed  [T],  If 
the  difference  between  any  corresponding  new 
and  old  values  exceeds  0.01  °C,  the  iteration  is 
stopped  and  the  old  temperatures  are  then 
replaced  with  the  newly-calculated  ones  and 
taken  to  be  the  required  temperatures  at  the 
section  concerned.  At  the  end  of  the  iteration, 
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the  progrmme  calculates  the  outlet  temperatures 
of  the  two  air  streams  and  also  the  efficiencies 
at  the  end  of  the  current  section.  It  then  sets 
the  air  temperatures  at  the  inlet  of  the  next 
section  equal  to  the  outlet  air  temperatures  of 
the  current  section. 

The  programme  then  proceeds  to  look  at  the 
next  section.  The  iterative  solution  procedure 
to  determine  the  new  mean  temperatures  of  this 
next  section  is  then  repeated.  By  this  repetitive 
and  iterative  process,  the  required  temperatures 
along  the  entire  length  of  the  solar  heater  could 
be  obtained,  section  by  section.  Normally,  the 
number  of  iterations  required  is  not  more  than 
two  or  three. 

5.  THEORETICAL  PREDICTION 

Values  of  absorptivity,  emissivity,  and  trans¬ 
missivity  for  the  various  surfaces  are  assumed 
values  taken  from  generally  accepted  figures 
and  not  measured  values.  From  preliminary 
results,  it  was  seen  that  small  temperature 
differences  of  the  order  of  1°C  or  less  could  be 
obtained  by  modifying  these  values.  This  meant 
that  the  prediction  was  not  heavily  dependent 
upon  these  values.  The  mathematical  model 
allows  for  different  values  of  heat  transfer 
coefficients  to  be  specified  for  the  upper  and 
lower  surfaces  of  the  walls  forming  the  air 
channels.  In  the  present  programme,  the  coeffi¬ 
cients  for  both  upper  and  lower  surfaces  of 
each  channel  were  assumed  equal,  hl  =  h2  and 
h3  =  h4. 

5.1 .  Effects  of  section  size 

As  stated  earlier,  the  theoretical  model 
assumes  that  for  a  short  collector,  the  wall 
temperatures  are  assumed  uniform  and  the  air 
temperatures  are  assumed  to  vary  linearly  along 
the  collector  in  the  direction  of  flow.  From 
observations,  these  assumptions  are  valid  for 
short  collectors.  However,  for  a  long  collector, 


the  temperature  rise  is  expected  to  be  asymp¬ 
totic.  In  order  to  determine  the  effect  of  section 
length,  three  different  sets  of  predicted  temper¬ 
ature  results  were  obtained  and  compared,  each 
with  a  different  value  of  section  length.  It  was 
observed  that  plate  temperatures  differed  by  less 
than  0.5°C  between  the  results  obtained  for 
0.5  m  and  1.0  m  sections.  Also,  the  air  temper¬ 
atures  and  efficiency  values  remained  nearly  the 
same  for  these  two  sets.  Hence,  for  the  present 
purpose  it  could  be  deemed  that  satisfactory 
results  could  be  obtained  by  using  aim  section 
collector. 

5.2.  Temperature  variation  along  flow  direction 

Predictions  were  obtained  for  the  four  types 

of  collectors  considered  by  assuming  air  mass 
flow  rates  equal  to  0.04  and  0.40  kg  s_1  for  the 
ambient  conditions  assumed  according  to 
Table  4.  For  discussion  purposes,  the  Reynolds 
Numbers  based  on  these  flow  rates  are  evaluated 
at  27°C  and  are  equal  to  4000  and  40,000. 

Wall  and  air  stream  temperatures  are  shown 
plotted  in  Figs  7-14  for  the  four  different  types 
of  collectors  at  Reynolds  Numbers  of  4000  and 
40,000.  These  correspond  to  flows  in  the  usually 
encountered  transition  and  turbulent  flow 
regimes.  It  can  be  seen  that,  as  expected,  as  air 
flow  rate  increases,  temperature  decreases.  Also, 
for  high  flow  rates  the  collector  operating  tem¬ 
perature  would  be  lower,  thus  resulting  in  higher 
efficiencies.  In  all  the  results,  it  was  observed 
that  the  absorber  plate  exhibited  the  highest 
temperature,  and  the  glass  the  lowest.  As  air 
flow  rate  increases,  the  temperature  differences 
between  the  air  stream  and  the  glass  bottom 
plate  decrease. 

5.3.  Effects  of  wind  convective  heat  transfer 
coefficient 

A  typical  value  of  radiation  heat  transfer 
coefficient  from  the  top  surface  of  a  collector  at 
100°C  and  ambient  at  27°C  from  eqn  (29)  is 


Table  4.  List  of  parameters  assumed  in  predictions 


Collector 

type 

Radiation 

intensity 

H  (W  m-2) 

Ambient 

temperature 

7a  (°C) 

Wind 

velocity 

V  (m  s-1) 

Air  flowrate 
stream  #1 
mi  (kgs-1) 

Air  flowrate 
stream  #2 
m2  (kg  s-1) 

1 

700 

27 

1.0 

0.04 

0.00 

1 

700 

27 

1.0 

0.40 

0.00 

2 

700 

27 

1.0 

0.04 

0.00 

2 

700 

27 

1.0 

0.40 

0.00 

3 

700 

27 

1.0 

0.00 

0.04 

3 

700 

27 

1.0 

0.00 

0.40 

4 

700 

27 

1.0 

0.04 

0.04 

4 

700 

27 

1.0 

0.40 

0.40 
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Distance  along  collector  (m) 

Fig.  15.  Effects  of  wind  heat  transfer  coefficient  for  Type  I  collector  at  Reynolds  Number  =  4000. 


01  23456789  10 

Distance  along  collector  (m) 

Fig.  16.  Effects  of  film  heat  transfer  coefficients  for  Type  I  collector  at  Reynolds  Number  =  4000. 
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about  10.5  Wm-2  K-1.  From  Table  1,  the  wind 
heat  transfer  coefficients  hw  from  the  different 
authors  varied  from  3.0  to  7.6  W  m-2  K-1  at  a 
typical  wind  speed  of 0.5 ms-1.  The  overall  top 
heat  loss  Ut  is,  thus,  appreciably  affected  by  the 
evaluation  of  hw.  The  effects  of  the  wind  heat 
transfer  coefficient  on  the  predicted  temperature 
results  are  shown  in  Fig.  15.  Two  sets  of  results 
for  the  Type  I  collector  at  a  Reynolds  Number 
equal  to  4000  are  compared  here.  The  first  set 
evaluates  hw  based  on  the  McAdams  (1954) 
correlation  whereas  the  second  set  is  based  on 
Watmuff  etal.  (1982)  correlation.  The  results 
show  that  by  decreasing  the  value  of  the  wind 
heat  transfer  coefficient,  and  thus  of  the  overall 
top  heat  loss  coefficient  value,  surface  and  air 
temperatures  are  appreciably  increased. 

Typical  rear  insulation  for  a  0.025  mm  thick 
polystyrene  sheet  (k  =  0.04  W  m_1  K_1)  is 
about  l,6Wm"2K_1,  Combined  with  the  wind 
coefficient,  the  overall  back  loss  coefficient  is 
not  substantial  and  would  not  have  any  appre¬ 
ciable  effect  on  the  predicted  results. 


5.4.  Effects  of  air-wall  convective  heat  transfer 
coefficients 

In  order  to  demonstrate  the  effects  of  the  film 
convective  heat  transfer  coefficients  on  the  over¬ 
all  results,  two  sets  of  predictions  were  obtained 
for  the  Type  I  collector  at  Reynolds  Number 
equal  to  4000.  The  correlations  listed  in  Table  3 
were  employed  in  one  set  and  in  the  other,  the 
resulting  value  of  the  heat  transfer  coefficients 
for  the  air  stream  was  multiplied  by  a  factor  of 
0.75.  The  predictions,  plotted  in  Fig.  16  shows 
that,  as  expected,  lower  coefficients  resulted  in 
a  decreased  cooling  effect  on  the  absorber  plate 
which  resulted  in  higher  plate  temperature  and 
lower  air  stream  temperature  and  thus  lower 
efficiency.  The  effect  on  the  bottom  plate  tem¬ 
perature  was  not  so  strong  as  for  the  others. 


6.  CONCLUSIONS 

This  paper  presents  a  theoretical  model  for 
predicting  the  thermal  performance  of  flat-plate 
solar  air  collectors  without  recourse  to  deriving 
complicated  equations.  The  mathematical  solu¬ 
tion  procedure  involved  a  matrix  inversion  of 
the  mean  temperature  vector  derived  from  the 
energy  equations.  Predicted  temperatures  for 
four  types  of  solar  air  heaters  were  presented. 


The  effects  of  wind  and  film  heat  transfer 
coefficients  on  the  prediction  were  discussed. 


NOMENCLATURE 

Ac  cross  section  of  flow  area  (m2) 

C  specific  heat  of  air  (1009.0  J  kg"1  K-1) 

Dh  equivalent  diameter,  4  x  flow  area/wetted  perimeter 
(m) 

Dhlh2  equivalent  diameter  of  channels  1  and  2  (m) 

/  friction  factor 

g  gravitational  constant  (9.81  m2  s  *) 
h,  2  3  4  forced  convection  heat  transfer  coefficients  (W  m  2 
K-1) 

hnc  natural  convection  heat  transfer  coefficient  (W  m  2 
K-1) 

hrb  bottom  plate — ambient  radiation  heat  transfer  coef¬ 
ficient  (W  m  2K-1) 

h„  top  cover — sky  radiation  heat  transfer  coefficient 
(Wm-2K-‘) 

/ir21r23  radiation  heat  transfer  coefficients  (W  m  2  K  1 ) 
hw  wind  convection  heat  transfer  coefficient  (W  m  2 

K-1) 

H  incident  solar  radiation  intensity  (900  W  m  2) 
k  thermal  conductivity  of  air  (W  m-1  K-1) 
kbi,b2  insulation  thermal  conductivity  (0.04  W  m-1  K-1) 
L  length  of  collector  (m) 

m12  mass  flow  rate  of  air  streams  1  and  2  (kg  :  *) 

P  wetted  perimeter  (m) 

Qu2  heat  transferred  to  air  streams  1  and  2  (W  m-2) 

S  factor  in  average  Nusselt  No.  evaluation 
Si  2  solar  radiation  absorbed  by  surfaces  1  and  2 
(Wm-2) 

ti  spacing  between  glass  and  absorber  surfaces 
(0.0254  m) 

r2  spacing  between  absorber  and  bottom  surfaces 
(0.0254  m) 

7^  2,3  temperatures  of  surfaces  1,  2  and  3  (K) 

T„  ambient  temperature  (K) 

Tfi,f2  mean  fluid  temperatures  (K) 

Ts  sky  temperature  (K) 

Tf  i  mean  fluid  temperature  at  inlet  to  section  (K) 

Tf  o  mean  fluid  temperature  at  outlet  of  section  (K) 

Ub  bottom  heat  loss  coefficient  (W m-2  K-1) 

Ut  top  loss  heat  coefficient  (W  m-2  K-1) 

V  wind  velocity  (1.0  ms-1) 

W  width  of  collector  (0.254  m) 
xbl  b2  insulation  thickness  (0.025  m) 

y  distance  along  collector  in  flow  direction  (m) 


Greek  symbols 


a 2  absorptivity  of  glass  (0.06) 
a2  absorptivity  of  black  absorber  upper  surface  (0.95) 
P  volumetric  coefficient  of  expansion  [2/(71  +  Ti) 
K-1] 

▲  T  temperature  difference  between  plates  [abs(T2  —  Ti) 
K] 

ej  emissivity  of  top  glass  surface  (0.90) 
e2  emissivity  of  black  absorber  upper  surface  (0.95) 
e3  emissivity  of  unpainted  absorber  lower  surface 
(0.25) 

e4  emissivity  of  unpainted  bottom  plate  surface  (0.25) 
H,  <"'2  efficiency  of  individual  air  stream  (%) 
n0  efficiency  of  combined  air  streams  (%) 

<t  Stefan — Boltzmann  constant  (5.67  x  10-8  W  m-2 


H  dynamic  viscosity  of  air  (kgm-1  s-1) 
<f>  tilt  angle  (7^°) 
p  density  of  air  (kgm-3) 
x  transmissivity  of  glass  (0.84) 
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